Negative Differential Resistivity and Positive Temperature Coefficient of Resistivity 
effect in the diffusion limited current of ferroelectric thin film capacitors 
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We present a model for the leakage current in ferroelectric thin- film capacitors which explains 
two of the observed phenomena that have escaped satisfactory explanation, i.e. the occurrence of 
either a plateau or negative differential resistivity at low voltages, and the observation of a Positive 
Temperature Coefficient of Resistivity (PTCR) effect in certain samples in the high-voltage regime. 
The leakage current is modelled by considering a diffusion-limited current process, which in the 
high- voltage regime recovers the diffusion-limited Schottky relationship of Simmons already shown 
to be applicable in these systems. 



For a number of years the problem of understand- 
ing leakage currents in ferroelectric capacitors (or high 
dielectric constant capacitors using ferroelectric materi- 
als just above their phase transtion, e.g. barium stron- 
tium titanate (BST)) has been of practical interest, 
and a large number of papers have been published on 
the subject, a review of which can be found in Scott's 
bookpj. Typically the currents observed appear to be ei- 
ther emission-limited or space-charge-limited. In this pa- 
per we restrict ourselves to the discussion of the emission- 
limited type currents. What this paper shows is that a 
metal-ferroelectric-metal system with space-charge due 
to oxygen vacancies can be considered as a diffusion- 
limited current system. The diffusion-limited current 
has two regimes: a low-voltage regime akin to the low- 
voltage regime discussed in the theory of semiconduc- 
tor punch-through diodes Q; and a high- voltage regime 
equivalent to the diffusion-limited Schottky relationship 
of Simmons 3] . 
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FIG. 1: Band diagram and representation of charges used to 
calculate field and potential distributions 

We consider the field distribution in a capacitor (Fig. 
1) which, as well as having a large dielectric constant 
£f and high uniform concentration of stationary oxy- 
gen vacancies Nd, has a remnant polarisation, P r . The 
thomas-fermi screening length in the electrode and the 



penetration of the metal states into the ferroelectric are 
accounted for by separating Q m , the charge on the metal 
from Q ss , the charge due to surface states (metal induced 
gap states) by an effective thickness <5 e //- This effective 
thickness reflects the fact that there is a finite screening 
length in the metal (the Thomas-Fermi screening length 
At/) and that metal-induced gap states penetrate a small 
distance into the insulator (Xmigs)- 



°eff = 



Xtf , Xmigs 
— H 



(1) 



Metal induced gap states have previously have been 
shown to be important in ferroelectric thin film capac- 
itors by Robertson and Chen 4] . The metal screening 
length term is typically the larger of the two, so S e f / is 

usually between about 0.5 and 1 ^ 

We now consider that the electric displacement in the 
film is given by 



D = e f s E±P r 



(2) 



The sum of the charges in the system must equal zero 
since the electric displacement is zero in both electrodes 
(far from the interface) . In the center of the film the field 
must equal the field that would arise if there were no sur- 
face states or space charge, which we take to be the linear 
drop = X-) that would occur if the sample were a per- 
fect insulator. This is quite different from the approach 
often taken, which is to treat a metal-semiconductor- 
metal system as two semi-infinite junctions, with the elec- 
tric field taken as zero in the regions of the film where 
there is no space charge, regardless of the applied field. 
Such an approach would be incompatible with the space- 
charge-free limit of a linear voltage drop for our insulat- 
ing system. (This approach is however perfectly appro- 
priate for a single infinite-metal/infinite-semiconductor 
junction). In addition, because our system is ferroelec- 
tric, when there is a spontaneous polarisation a finite 
potential drop across the system occurs because of the 
imperfect screening of the depolarisation field, even in 
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the absence of an applied field. Using these arguments 
one finds that the field in the space-charge-free interior 
of the film is 



En — — ■ 



1 V±P r 2S eff 
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Our expressions for the charge on the metal electrodes 
are then 



Qm 1 = -Q ssl - qN D wi - — '-T- (4) 
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Qm 2 = -Q SS 2 - qN D w 2 + 



2<5, 
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(5) 
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The amount of charge due to surface states is inde- 
pendent of voltage and is found by the condition that 
the fermi level at the interface of the ferroelectric is 
S{(j> m )) + (1 — S)(x + (E c - (/))) below the vacuum level 



where S = 
this gives 



l+q 2 5e ff D„ 



and 4> is the charge neutrality level; 



Q ssl = (l-S)(-qN DWl + — ((<j> m - X )-(E g -<j>))) (6) 

°eff 

This result shows that the charge due to surface states 
is independent of voltage, and under the assumptions 
that the two electrodes are the same and that the film is 
uniform, Q ssl = Q ss2 . 

Finally, the width of the depletion layers is found using 
the condition that the potential at x = wi must be equal 
to the height of the donor trap level E t (which in samples 
with low intrinsic conductivity and high dopant levels will 
be the semiconductor Fermi level). It is found that the 
depletion regions are invariant with applied voltage, 



Wl + SfSaSSeff = [((efEnSSeff) 2 + 



2e/£o 
qN D 



((0m - X) 



- (E c - E t ) - (1 - 5)((0 m - X ) - (E c - 0))]* (7) 

and that for a symmetrical system w 2 = w\ . If a film is 
fully depleted [i.e., if the calculated values of w\ and w 2 
are greater than ^ (which means that the system fermi 
level never reaches the trap level)], the following equa- 
tions may be used with w\ = w 2 = ^. 

The finding that the depletion widths do not change 
with field is at first quite surprising, the typical idea be- 
ing that in punch-through diodes, depletion widths will 
change with field and that at a certain applied bias the 
reverse-biased depletion width will "punch-through" to 
the forward-biased Schottky barrier. But in fact what 
we show here is that it not the depletion width itself that 



changes size, but rather the physical location within the 
system of the field inversion point. It is when this field in- 
version point coincides with the forward-biased Schottky 
barrier that punch-through behaviour occurs. 

The potential distributions for < x < w\, Wi < x < 
d — w 2 and d ~ w 2 < x < d respectively are V(x) = 



QmlOeff - E 1 X(W-L - -), 

Ef£ £/£o 2 
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QmlOeff ~ E h W 1 - [X -(d-W 2 )j 
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which give the field distributions in the same regions 



as 



E(x)=E -^( Wl -x), 



£/£ 

E(x) 
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E{x)=En + ^{x-(d~w 2 ) 
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(9) 



An important result of this system of equations is that 
below a certain critical voltage there is a point 



xo = - ^rr-En + d-w 2 
qN D 



(10) 



in the film at which the field is zero (the field inversion 
point referred to above). We can calculate the current 
flowing in the device by determining the point in the film 
at which the field is zero. Here the current is entirely due 
to diffusion and is given by 



ox 



(11) 



To evaluate §P- we look at the number of conduction 

ox 

electrons in the region d — w 2 < x < d, 



n(x) = N C exp(--^((£ c - E t ) + ^-(x (d-w 2 )f)) 

(12) 

Taking the derivative and using the solution for xn 
from above we find the current to be 



J = »N C E exp(-^((S c - Ek) + ^(Eo) 2 )) (13) 
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An equivalent expression for the current, which is also 
applicable for fully depleted films (in contrast to the ear- 
lier expression, which applies to only partially depleted 
films) is 



J = fiNcE exp(— — ( 



qNrjW 2 



e f e 



2s f e 2qN D 



(14) 



For a fully depleted film one simply substitutes w = i 
into the equation above. These relationships are applica- 



ble until xq = d which is when 



= qN D W2 [or 



2S e ffEf£o+d 

in other words when the field at the cathode is zero] . This 
is equivalent to punch- through, but a punch-through of 
the field inversion point, rather than the edge of the de- 
pletion width; and it marks the transition to a new con- 
duction regime. The current at which the transition oc- 
curs is the same as Frank and Simmons obtainjfj, i.e., 



qN D w 2 , <pb s n 
J T = -fiN c exp (- — ) 15) 

£/£0 Kl 

The voltage at which the transition occurs, 



V = qN D w 2 (26 eff + ) (16) 

£/£ 

is temperature- and thickness-dependent, the temper- 
ature dependence is somewhat complicated because of 
the temperature dependence of the dielectric constant 
(and hence also the depletion width). When the applied 
voltage is increased beyond this voltage, the image force 
pushes the potential maximum back into the film, low- 
ering the effective barrier. If we assume that the new 
potential maximum is quite close to the interface, we can 
use the field at the interface to calculate the Schottky 
barrier lowering as 



Q 



47T£ s £o 



(-Eq- 



qN D w 2 , 

£/£ 



(17) 



The dielectric constant £ s used here is the optical di- 
electric constant of the material and is used because elec- 
trons passing through the Schottky barrier do so suffi- 
ciently quickly that they do not polarise the lattice 6] . 

In this regime the current is 



J = - t ,NcE cxp(-j^(fo - A J) 



(18) 



We note that this is the form of the Richardson- 
Schottky equation applicable in materials with very short 
mean free paths and has been previously derived by Sim- 
mons H . This form of the Richardson-Schottky equation 
has been already been noted as the appropriate one for 



ferroelectric thin films and has been used successfully to 
describe leakage current in BST thin films by Zafar[7J 

In a previous study on Au-BST-SrRuC>3 capacitors^ 
we incorrectly assumed that the high- voltage regime cor- 
responded to tunnelling through the space-charge region 
in the electrode because of the correspondence between 
the threshold for the high- voltage regime and the applied 
voltage at which the potential drop across the electrode 
equalled the barrier height; i.e., it was noted that 



V T = 



2S, 



eff 



"( 2( W/ 
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(19) 



We note however that the present model predicts the 
same thickness dependence of the threshold voltage: V oc 
(25 e ff + jj^-)- Low temperature measurements on the 
70 nm film from this study show that the high-voltage 
regime is Schottky-like both at room temperature and at 
70 K (Fig. 2). 
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FIG. 2: Leakage current data from 70nm thick Au-BST- 
SrRuO?, film at room temperature and at T—70 K 

The low-field expression derived above (Eqs. 15 and 
16) can account for two kinds of low-field behaviour ob- 
served in ferroelectric thin films: either a sharp increase 
followed by a heretofore unexplained flat plateau (which 
is the behavior of the equation when the donor concen- 
tration is very high) or a negative differential resistivity 
(observed for lower donor concentrations). 
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FIG. 3: (a) Scott et al. r 9J: Leakage current of BST thin film 
showing fiat plateau followed by Schottky regime, (b) Watan- 
abe et al. \l Leakage current of SBT film showing negative 
differential resistivity followed by Schottky regime 
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An example of the first kind of behavior (Fig. 3(a)) 
is found in the data of Scott et al 9] on BST thin films; 
whereas a good examples of the second (Fig. 3(b)) can 
be found in the SBT samples of Watanabe et alllCi or 
in PZT in the data of Scott et al.[ll| or Chen et al|l2|. 
Previous attempts to explain these effects have invoked 
the filling and emptying of trap states. 

On taking the derivative with respect to voltage of our 
low-field expression, we find that 



dJ_ 

dV 



riN c 



(1 - (£o) 



2 £ f £ 



kTN 



D 



exp(--^(</>6 - 



qN D w 2 
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(20) 



The exponential term decreases with voltage towards 
the value exp(— ^) but is always positive. If the other 
voltage-dependent term remains positive, the low-voltage 
leakage current characteristic will be a rapid increase 
which plateaus off into an ohmic regime until such volt- 
age that the barrier begins to be lowered by the Schottky 
effect. However, if the voltage is greater than 



V > (2S ef je f £ + d) y/kTej^No (21) 

the differential resistivity will become negative. We 
should not forget, however, the condition that this con- 
duction mechanism only operates when 



V < 2qN D w 2 (25 ef f + ) 

£f£ 



(22) 



and so negative differential resistivity is observed only 



if 



/ e/epfcT 
Aq 2 N D wl 



< 1 



(23) 



Within the Schottky injection regime one of the most 
interesting things to note is the role that the dielectric 
constant plays in determining the temperature depen- 
dence of the current. We can show that when the sample 
is below a ferroelectric phase transition and the dielec- 
tric constant is increasing rapidly with temperature a 
regime of Positive Temperature Coefficient of Resistivity 
(PTCR) may be observed. For this purpose we simplify 
the Schottky regime current expression by taking 
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(-Eo) 



(24) 



and taking the derivative with respect to temperature 
we find that 



— = J(--^-\-( 
dT v kT T 



6 b -A0 fc )-A0 b ( 



d£ f ( 

dT { 25 eff £f£o + d 



J eff 



))]) 



The term ^ ( 



(25) 

bjj — Atpb) is always positive. For a PTCR 
effect to be observed above a certain applied voltage 
d< dT° ( 26 /jejeo+d ) mus ^ be large and positive, and hence 
the effect should be seen only below Tc- Whether or 
not PTCR behaviour is observed depends on the barrier 
height and the screening length in the metal, with PTCR 
behaviour more likely to occur for capacitors with lower 
barrier heights and longer screening lengths. The PTCR 
effect predicted here does not involve grain boundaries, 
as distinct from the PTCR effect in BaTiC>3 ceramics (for 
a review see Huybrechts et al. jj^) which as shown by 
Sinclair and West has both bulk and grain bound- 
ary contributions. Hwang 15] has observed PTCR effects 
in films with Ir electrodes, whereas similar films with Pt 
electrodes did not display PTCR behaviour. This is con- 
sistent with our model because Ir has a significantly lower 
workfunction than Pt (4.23 eV compared to 5.3 eV [U|). 

Thus with a single model we have explained many un- 
usual leakage current features in ferroelectric capacitors 
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